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Abstract 

Using a geometry-based fundamental measure density functional theory, we calculate bulk fluid 
phase diagrams of colloidal mixtures of vanishingly thin hard circular platelets and hard spheres. 
We find isotropic-nematic phase separation with strong broadening of the biphasic region upon 
increasing the pressure. In mixtures with large size ratio of platelet and sphere diameter, there is 
also demixing between two nematic phases with differing platelet concentrations. We formulate a 
fundamental-measure density functional for mixtures of colloidal platelets and freely overlapping 
spheres, which represent ideal polymers, and use it to obtain phase diagrams. We find that for low 
platelet-polymer size ratio in addition to isotropic-nematic and nematic-nematic phase coexistence, 
platelet-polymer mixtures also display isotropic-isotropic demixing. In contrast, we do not find 
isotropic-isotropic demixing in hard core platelet-sphere mixtures for the size ratios considered. 
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I. INTRODUCTION 



Understanding the equilibrium properties of colloidal systems and the relationship be- 
tween the microscopic properties, such as particle shapes and sizes, and the macroscopic 
properties of a dispersion is essential for the task of creating new materials with desired 
characteristics. Here anisotropic colloids are relevant due to the possible formation of liquid 
crystalline phases. Among them platelet-like colloids are ubiquitous in nature; gibbsite or 
certain clays are specific examples. 

Despite their simplicity, hard-core models are very suitable candidates to investigate the 
phase behaviour of colloidal systems. Since Onsager's pioneering work on the isotropic- 
nematic (IN) phase transition in a fluid of rods a wealth of studies has been carried 
out in order to elucidate the equilibrium properties of anisotropic hard particles, including 
the fluid of circular platelets with vanishing thickness using Onsager-like theories 3| and 



fundamental measure theory (FMT) Q,H|; see Ref. 



for a recent review of the literature 



for platelets and Refs. |7H9l for reviews of FMT. 

During the last few years, binary colloidal mixtures have received considerable attention. 
The addition of a second component induces profound changes in the phase diagram with new 
phenomena arising such as demixing, reentrant phase boundaries, as well as the emergence 
of phase transitions due to the depletion mechanism. Studies of binary mixtures that involve 



platelets were conducted for mixtures of thick and thin platelets 



with different diameters 



12l |. and mixtures of platelets and rods |13L \l4 \ 



Ml 



Jj, binary platelets 



In this paper we use density functional theory (DFT) [l5|, and in particular a geometry- 
based fundamental measure free energy functional Q|, in order to investigate the bulk fluid 
properties in mixtures of infinitely thin platelets and hard spheres. Such mixtures have 
been previously considered by Harnau and Dietrich using DFT {l6| and by Oversteegen and 
Lekkerkerker by means of a free- volume scaled-particle approach 17|. The authors restricted 
their analyses to isotropic states of platelets. Both studies found coexistence between an 
isotropic fluid rich in platelets and a solid phase rich in spheres, as well as the stability of a 
phase transition between two isotropic phases with differing compositions for small values of 
the platelet-sphere size ratio. Both studies neglect the excluded volume interaction between 
the platelets, which can certainly be a good approximation if the density of platelets is very 
low. We will compare to these findings below. We compare also to results from an extended 
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Onsager theory with Parsons-Lee rescahng 



19| which includes two-body correlations 



only. In contrast, the geometry-based density functional that we use here includes explicitly 
non-local higher-body correlations. The non-local structure of the free energy functional 
is specifically tailored for the present shapes. For general convex bodies much recent effort 
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22| was devoted to improve Rosenfeld's original FMT for non-spherical shapes 
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Here we also study mixtures of platelets and non-adsorbing ideal model polymers 



28 



26|. 



I.e., 



Polymers are modelled following the ideas of Asakura and Oosawa |27|, and Vrij 
they are treated as fully interpenetrable spheres that cannot overlap with the colloidal 
platelets. This is to be regarded only as a minimal model, in particular for the platelet- 
polymer interaction, as both softness and penetrability of a platelet into a polymer coil are 



neglected. We develop a fundamental measure density functiona 
density functional for hard-core mixtures in polymer density 29 . 



theory by linearizing the 



3Q|. 



Experimental mixtures of sterically stabilized gibbsite platelets and silica spheres have 
been studied very recently jsil- The authors found phase separation into platelet-rich and 
sphere-rich phases, driven by the repulsive interactions between particles of different species. 
Columnar ordering has been observed in charged mixtures of gibbsite platelets and silica 



spheres 



32 



33] . Furthermore the crystallization of silica spheres in presence of gibbsite 



platelets (silica-coated) was analyzed in Ref. 34| . 

The paper is organized as follows. In Sec. [Illwe define the model and outline the DFT. The 
results are described in section lllli mixtures of colloidal platelets and spheres are analyzed 
in Sec. IIII Al in Sec. IIII BI we present the results for mixtures of colloidal platelets and ideal 
polymers, along with a comparison between FMT theory and an extended Onsager theory 
which is given in IIII CI Suggestions for future research lines and conclusions are presented 
in Sec. IVl 



II. MODEL AND THEORY 



We consider a binary mixture of Ns hard spheres with radius Rs and N^) hard circular 
platelets with vanishing thickness and radius Rd, see Fig. [1] for a schematic illustration. 
Here and throughout species-dependent quantities are indicated by subscripts S (spheres) 
and D (disks). The pair interaction potential between any two particles is infinite if the 
particles overlap, and vanishes otherwise. 
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FIG. 1: (Online version in colour.) Schematic representation of the hard core platelet-sphere 
mixture. The model is a binary mixture of hard spheres of radius Rs (right) and hard circular 
platelets of radius Rd with vanishing thickness (left). The platelet orientation is described by a 
unit vector Cl normal to the platelet surface. 

A. Density functional theory for colloidal platelet-sphere mixtures 

As is common, we split the Helmholtz free energy functional F into an ideal-gas part, 
Fid, and an excess contribution, Fg^c- Hence 

F = Fid + Fexc- (1) 

The ideal-gas free energy functional is exactly given by 

/3Fid= [ dh [ dnpD{r,mn{pD{r,n)VD)-l] 

(2) 

+ J d'rps{f)[Hps{f^Vs)-l], 

where /3 = l/ksT is the inverse thermal energy, ks is the Boltzmann constant, T is absolute 
temperature, and is the (irrelevant) thermal volume of species i. The spatial integral is 
over the system volume, V, and the angular integral is over the unit sphere; here f2 is a 
unit vector normal to the platelet surface, cf. Fig. [1] The one-body distribution function 
of platelets is denoted by p£,{f,Cl), and that of spheres by P5(r), the latter being trivially 
independent of orientation due to the rotational symmetry of the spheres. In general both 
distribution functions depend on position r. In the application to bulk fluid phase behaviour 
below we consider only homogeneous states, which are independent of r. Without loss of 
generality we introduce an orientational distribution function for the platelets, \I^d(^, via 
pij{f,Cl) = pj:)(r*)\l/£)(r, (i). Here Pd(^) is the number density of platelets, which gives the 
differential number of platelets with any orientation in the differential volume element at r. 



The excess part of the free energy functional, Fexc-, that accounts for the excluded volume 
interactions between all particles, is approximated by a geometry-based density functional, 
which has the structure 

/3F,.e = j rfV$(K}), (3) 

where the reduced free energy density, $, is a function of a set of weighted densities, {n^}, 
where i = S,D labels the species, and u the type of weight function. The free energy 
density can be grouped into three different contributions, $ = $5 + $/5 + ^sd- Here $5 
describes sphere-sphere interactions only. Both and $5/5 describe coupling of spheres 
and platelets, with containing the terms that generate the FMT functional for pure 
platelets |4| in the limit of vanishing density of spheres. In detail, the reduced excess free 
energy for the interaction between spheres, $5, is taken to be of the original Rosenfeld form 



23|, 



$, = -ngln(l-nf) + "^"^ "-^ 

" (4) 

247r(l - nf )2 ' 

Here and in the following two equations we drop the spatial dependence on r of the weighted 
densities for the sake of notational convenience. The interactions between platelets are 
described by 



'-"^ .... (5) 
+ JdilJdU 24vr(l-nf)2 
The first term in Eq. represents the interaction between two platelets. It reduces to the 
exact second order contribution of the virial series for the free energy of platelets if there 
are no spheres in the mixture [ps = 0). The second term corresponds to the interaction 
between three platelets (see Ref. % ll2| for details). The factors proportionals to (1 — nf ) 
in the denominators account for the volume occupied by the spheres in the system. 

The further contribution to the reduced excess free energy due to the interaction between 



5 



platelets and spheres is 

^SD = - j dnn^{n) ln(l - 



3 J 

1-nl 
87r(l-nf)2 



(6) 



87r(l-nf)2 

All weighted densities are obtained by convolving the density with specific weight func- 
tions. For hard spheres the scalar weighted densities are: 

n^A^ = w'Af)*psif), z/ = 0,l,2,3 (7) 

and the vectorial weight densities are given by: 

<(f) = w'^if) * psif), V = vl, v2, (8) 

where * denotes the three-dimensional convolution h{r) * g{r) = ( d^xh{x)g{x — r). Here 
and are the Rosenfeld weight functions for hard spheres |23| : 

wl{r) = Q{Rs-\r\), 

w^{r) = 5{Rs-\r\), (9) 

where 6(-) is the Heaviside step function and 6{-) is Dirac's delta distribution. Further 
weight functions for hard spheres are wf{f) = (r)/(47ri?5), Wo(r) = (r ) / (4-7ri?|) , and 
(r-) = w^^{f)/{A'KRs). 

The sphere-platelet coupling is described by the weighted densities 

nf^(f,fi)=w;f^(f,fi)*p5(r), 

n) = w^i'ir, Q) * p5(r), (10) 

where the weight functions are 



wnf,n) = -Ji-if-n/RsrwUr), 

<F{r. fi) = /~^"-^^" ^f(f), (11) 

TlRs 



,SDDi 



71 
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Finally, the weighted densities for platelets are 



n 



DDi 



(f, n') = / dQw^^if, poif, ft), 



and their corresponding weight functions are 



1 



27tR 



D 



-SiRn-\n)5if-n), 



w^{f,n) = i-siRD-\r\)sif-n), 



w^{f, n) = 2<S){Rd - \r\)6{r- fi), 



w 



DD, 
1 ' 



f, fi, VL') = 



R 



D 



TV 



(12) 



(13) 



The density functional described here is a special case of the functional for ternary 
platelet-sphere-needle mixtures j4|. It includes contributions up to third-order in platelet 
density. It reduces to the original Rosenfeld functional 23] in the case of a monodisperse 
fluid of hard spheres, and gives results that are in good agreement with Monte Car 
lation data when describing the monodisperse fluid of vanishingly thin platelets j4. 
Hence we are confident that it can be applied successfully to the study of binary mixtures 
of hard spheres and platelets. 
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B. Density functional theory for mixtures of colloidal platelets and ideal polymers 

We use a model |26| similar in spirit to the colloid-polymer model of Asakura and Oo- 



sawa |27[, and Vrij [28|, to study a mixture of colloidal platelets, which are as above taken as 
infinitely thin circular disks, and ideal polymers. The polymers are modeled by fully inter- 
penetrable spheres of radius Rp. These model polymers cannot overlap with the platelets. 
We obtain an approximation for the free energy of this system by linearizing the excess part 
of the free energy for the corresponding hard core system. This strategy follows the deriva- 
tion of the FMT for the Asakura-Oosawa model from a corresponding binary hard sphere 



mixture |29l . |30j. Hence we apply the linearization to the hard core functional of Sec. IIII Al 
in order to obtain a free energy functional for the polymer-platelet mixture, F^J;^ , via 

Fex^ [pD, PP] = Fexc [PD, PS = O] 



(14) 

ppir), 

ps=o 



where pp{r) is the one-body polymer density distribution and F^xc is the free energy func- 
tional for the hard core mixture, given in Eq. F^J^ can be expressed in terms of a 
reduced free energy density 

PKx^'c = I rfV$pz,({<}), (15) 

where the reduced free energy, $pd, is a simple expression linear in polymer density and 
with up to cubic contributions in platelet density: 



di / ^n-"'°°'"-">?'»'"?'"'' (16) 



Stt 

+ / dnn^^{n)n^{n) (i + nf) 



J 24:71 V 3; 

We have suppressed the spatial dependence of the weighted densities in the above expres- 
sions. The weighted densities that depend on the polymer density (denoted by the letter 
P in the superscript) are those for hard spheres ([7]) and ([8]), but replacing the hard sphere 
colloid density, psir), by the polymer density pp{f), and replacing the colloid radius, Rq, 
by the polymer radius Rp. 



C. Spatially homogeneous fluids 

In the following we restrict ourselves to considering spatially homogeneous fluid states. 
We use the sphere-platelet model, the expressions for the polymer-platelet model being 
analogous. The density distributions are independent of the spatial coordinates and 

psir) = ps = const, 

(17) 

PD{f,Cl) = p£)\l/ £)((]), where po = const 
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The total number of spheres, Ns, and that of platelets, Nd, are obtained by integrating the 
number densities over the system volume 



Ns = / d-'rps, 

% (18) 

Nd = / d^rpD. 



Eq. (fT8|) implies that the orientational distribution function of platelets is normalized ac- 
cording to 

j dn^Di^) = 1, (19) 

which is a different convention compared to that of Ref. As we do not expect biaxial 
arrangements of the platelets to occur in bulk, we take the orientational distribution function 
to depend only on the polar angle 6 with respect to the nematic director, hence "^{Cl) = "^{0). 
We use a parameter. A, and prescribe the functional form of the orientational distribution 
function as 

^ = exp[AP,(cosg)] 
^ /dl]exp[AP2(cos^)] 
satisfying by construction the normalization condition (fT9|) . Here P2(") is the second Leg- 
endre polynomial, and A determines the degree of orientational order of the platelets. For 
isotropic states "^dI^) = l/47r is obtained for A = 0. We characterize the orientational order 
of the platelets by Sd, the standard nematic order parameter: 

Sd = j dh<i!D{0)P 2{cos 9). (21) 

We have tested the quality of the parameterization ( l20l) by computing the IN phase 
transition of the monodisperse fluid of platelets, and comparing the coexisting densities 
and order parameter to those obtained from free minimization Q]. The differences for 
the coexistence densities between both methods are of the order of 1%, see Table [T] for a 
comparison. The value of the nematic order parameter at coexistence, S^, deviates more 
strongly from the result of free minimization. However, as its magnitude is unusually small 
for the pure platelet system, this is a very sensitive quantity, see the discussion in Ref. |l2l |. 
For completeness. Table |T] also gives the results from Onsager theory and from Monte Carlo 
simulations. 

In spatially homogeneous fluids, the weighted densities are obtained by integrating the 
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phRo 


PdK 


C'i V 


FMT, free minimization [121 


0.419 


0.469 


0.533 


FMT, parametrization ([20]) 


0.421 


0.462 


0.489 


Monte Carlo [36| 


0.460 


0.498 


0.45 - 0.55 


Onsager theory 


0.672 


0.858 


0.801 



TABLE I: Coexistence values for the IN phase transition of a monodisperse fluid of hard circular 
platelets with vanishing thickness. The values in the first row are from free minimization without 
an a priori form of the orientational distribution function |l2|. T he second row shows the coexisting 



36| are presented in the third row. In 



values using the parameterization (|20|) . Monte Carlo data 
the fourth row the results from Onsager theory are shown; these are identical to Parsons-Lee theory 
(Sec lIIP]) . as the packing fraction vanishes in the pure system of platelets. Here, and are the 
coexistence densities at the IN transition, and is the nematic order parameter in the coexisting 
nematic state. 



weight functions over the spatial coordinates. For iiard spiieres: 

'^o = Ps, 



nf = Rsps, 

nl = AnRlps, (22) 
<i = <2 = 0. 



The mixed weighted densities that couple spheres and platelets are: 



(23) 



nl^^{9,9') = Ps dz dcf)' I (24) 



-Rs Jo \ 64 (Vs2 - t2 + t 

arcsin -) otherwise 

where t = z sin 9 sin 9' sin ch' , and 



s = JrI - z'^J sin^ 9 cos^ 9' + cos^ 9 sin^ 9' - 2 sin 9 cos 9 sin 9' cos 9' cos 0'. (25) 
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Due to our normalization ( IT9|) . Eq. above and Eq. (105) of Ref. (j| differ by a factor 
(47r)^. The integrals over the azimuthal angles (p and (p' have been performed in order to 
arrive at (12^ . as only the polar angle with respect to the nematic director is relevant for 
uniaxial configurations. For platelets the bulk weight functions reduce to 

4 (26) 



nf^m = 0. 



Further contributions to the free energy due to the weighted densities nf ^ and n^^^ are 



the same as given in Eq. (13) and (27) of Ref. |l2j to which we refer the reader directly. 

The angular integrals were calculated numerically using Gaussian quadrature. We used 
2000 roots for the integrals over the azimuthal angle and 25 roots for the integrals in the 
polar angle. The orientational distribution function does not depend on the azimuthal angle, 
and hence the integrals over can be computed and stored for subsequent use. We checked 
the accuracy of the numerical procedure by computing selected coexisting points with 100 
roots for the integrals over the polar angle. The differences were less than 1%. 

D. Parsons-Lee theory for colloidal platelet-sphere mixtures 

Parsons-Lee theory 



19] constitutes an approximation where the second virial contri- 
bution to the excess free energy is weighted with a prefactor in order to approximate the 
effects of higher virial contributions. According to Parsons-Lee theory the excess free energy 



for a spatially homogeneous binary mixture of hard bodies can be written as (see e.g. 37|): 
^ = m YIp^P^ j'^^j d^'^m'^,{n')v%{^l,^l'), (27) 

where v'^^^{VL, VL') is the excluded volume between two particles of species i and j with orien- 
tations given by f2 and fi', respectively. The prefactor 0(?7) depends on the chosen reference 
system. Here rj is the total packing fraction of the mixture rj = YliPi'^i^ where Vi is the 
particle volume of species i. Our reference system is a pure fluid of hard spheres, and for 
the present model rj = Vsps, as the platelet volume vanishes. Using the Carnahan-Starling 
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equation of state, the prefactor in the excess free energy is 

The orientational distribution function for spheres is = l/4vr, and for platelets is 

given in Eq. (pOl) . The excluded volume between two spheres is 

v!^{^,^') = ^^R% (29) 

between platelets and spheres 

^;f/((i, h') = 27iRsRl + tt^RdRI + l^Rl (30) 

and, finally, between two platelets 

vg''{n,n')=AnRlsmj, (31) 

where 7 is the angle between Cl and Q'. For pure platelets ?7 = 0, and hence (128|) reduces to 
(f){ri) = 1/2, which renders fl27|) identical to the (Onsager) second-virial functional. 

E. Coexistence conditions 

In what follows we denote the composition of the mixture, x, by the molar fraction of 
hard spheres, i.e., x = ps/p, where p = ps + Po the total density. The molar fraction 
of platelets is then simply p^/ p = I — x. The equilibrium properties of the mixture are 
obtained by minimizing the Gibbs free energy per particle g = F/Nt + p/p at constant 
composition x, pressure p, and temperature T. Here F is the total Helmholtz free energy, 
and Nt = Ns + N^, is the total number of particles. We use a standard Newton-Raphson 
method to minimize (7 as a function of p and A. Binodal lines are located by a common- 
tangent construction on g{x), which is equivalent to the equality of chemical potentials of 
both species at the coexisting values of x and p. Thermal and mechanical equilibrium are 
satisfied in advance by fixing both temperature and pressure. 

Correspondingly, for the platelet-polymer mixture, x = pp/p, and p = pp + po- The 
coexistence conditions are those given above with colloidal spheres quantities exchanged by 
the respective polymer quantity. 
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III. RESULTS 



The pure fluid of hard platelets with vanishing thickness undergoes a first order IN phase 
transition as a function of density. Table [T] summarizes the values for the coexisting densities 



and for the nematic order parameter according to t 
agreement with Monte Carlo simulation data 38 



le present FMT. The results are in good 



39|. Columnar and solid phases appear 



only in the limit of infinite density in a pure system of platelets with vanishing thickness 
|36|. In a monodisperse fluid of hard spheres there is a well-known first-order fluid-solid 
phase transition by increasing the density. We have not studied spatially ordered phases, 
but give some comments about the stability of spatially homogeneous phases with respect 
to hard sphere crystallization in Sec. IIV[ 



A. Mixtures of colloidal platelets and spheres 

In Fig. |2]we present the phase behaviour of mixtures of platelets and spheres for three dif- 
ferent mixtures with size ratios Ro/Rs = 2 (first column), 5 (second column), and 10 (third 
column). We comment on the behaviour for smaller size ratios at the end of the section. In 
Fig. [2] (al) we show the phase diagram in the plane of reduced pressure, fSpR^, and molar 
fraction of spheres, x, for a mixture with size ratio R^/Rs = 2. The IN phase transition in 
the pure fluid of platelets (x = 0) moves to higher pressures and the fractionation becomes 
stronger when spheres are added to the system. Upon increasing the pressure, the nematic 
branch of the binodal initially moves to higher compositions, but then bends back on itself 
approaching x = 0. The nematic phase is stable only when the molar fraction for spheres is 
very low. 

In Fig. [2] (a2) we plot the phase diagram for the same mixture in the plane of reduced 
density of disks, poR^, and packing fraction of spheres, psVs, where vs is the sphere volume. 
The IN biphasic region strongly broadens when the packing fraction of spheres in the 
mixture is increased. The orientational order parameter of platelets, Sd, along the nematic 
side of the binodal is shown in panel (a3) of Fig. |2l We also display the total orientational 
order parameter of the mixture i.e., the nematic order parameter of platelets weighted by 
the composition of platelets: 

Stot = Pd/pSd = (1 - x)Sd- (32) 
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FIG. 2: (Online version in colour.) Phase diagrams of binary mixtures of spheres and platelets with 
three different values of the platelet-sphere size ratio: Ro/Rs = 2 (first column), 5 (second column), 
and 10 (third column). For each mixture we show the reduced pressure vs composition of 

spheres x (first row), the reduced density of platelets poRf) as a function of the packing fraction for 
spheres psvs (second row), and the nematic order parameters along the nematic side of the binodal 
(third row). The shadowed areas indicate two-phase regions. Empty circles represent critical points. 
Dotted lines represent triple points. The dashed lines in the second row are selected tie lines that 
connect two coexisting points on the binodal. The inset in panel (b2) is a zoom of the region close 
to the IN phase transition of the pure fluid of platelets. The solid black (dashed red) lines in the 
third row represent the platelets' (total) nematiji^prder parameter, So (St)- The insets in the third 
row are cartons showing the relative size between species. 



As the composition of spheres at the IN phase transition is low, both orientational order 
parameters are similar to each other, increasing along the nematic side of the binodal. 

The second column of Fig. [2] shows the results for mixtures with R^/Rs = 5, i.e., for 
increased relative size of the platelets. The topology of the phase diagrams in both represen- 
tations is the same than in the above mixture. The main difference to the case of R^i/Rs = 2 
is the strongly increased range of compositions for which the nematic phase is stable. For 
example, in Fig. [2] (a2) one can observe that, for a range of pressures, it is possible to find a 
stable nematic phase even if the composition of the mixture is x ^ 0.6. The mechanism that 
induces the IN phase transition is different at low and at high pressures. At low pressures 
the transition is mainly due to the gain in configurational entropy of the platelets in the ne- 
matic phase, as the excluded volume between two platelets is minimal if they are parallel to 
each other. In this regime there is low partitioning between the isotropic and nematic phases. 
At high pressures, however, the transition is driven by the (unfavourable) excluded volume 
interactions between platelets and spheres. As a result there is strong demixing between a 
nematic phase rich in platelets and an isotropic phase mostly composed of spheres. 

We focus next on the behaviour of p/j along the nematic branch of the binodal. For low 
size ratios it monotonically increases as spheres are added to the mixture; but for Rs ^ 2, 
Pd slightly decreases when the packing fraction of spheres is low (see the inset in panel (b2) 
of Fig. |2]). This behaviour appears to be counterintuitive. However a reduction of the density 
of platelets at which the nematic phase is stable in the mixture has recently been observed 
in sedimentation experiments on sterically stabilized mixtures of gibbsite platelets and silica 
spheres jsil of size ratio ^ 6. However, the authors found a much stronger effect than that 
obtained in our theory. The quantitative difference could be due to the influence of gravity 



on the bulk phase diagram 



40|. The polydispersity of the components, the finite thickness 



of the platelets, and further interactions between the particles (not purely entropic as we are 
considering here) in the experimental setup could also be relevant to explain the quantitative 
differences. 

The last column in Fig. |2] shows results for the phase behaviour of a mixture with size 
ratio Ro/Rs = 10. A prominent alteration of the topology of the phase diagram is that 
in addition to the IN phase transition, there is also demixing between two nematic phases 
with differing compositions. The nematic-nematic phase coexistence region is bounded by a 
lower critical point and merges with the IN phase transition at an isotropic-nematic-nematic 
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triple point. At pressures above the triple point pressure there is strong demixing between 
an isotropic phase rich in spheres and a nematic phase mostly composed of platelets. The 
nematic-nematic demixing can be viewed as being driven by a depletion attraction induced 
between the platelets by the small spheres. 

The behaviour of the orientational order parameters along the binodal is shown in the 
third row of Fig. [2l For small size ratios, panel (cl), both the order parameter of platelets, 
Sd, and the total order parameter, Stot, increase along the nematic branch of the binodal. 
However, in mixtures where the size ratio is high (panels (c2) and (c3) of Fig. |2]) Sd mono- 
tonically increases along the binodal, but the total order parameter, Stot, first decreases to a 
minimum and then increases. This qualitative difference may be observable experimentally. 

The phase diagram of mixtures with size ratios Ru/Rs < 2 (not shown) are qualitatively 
similar to the results for Ro/Rs = 2. In particular, the topology remains unchanged and 
only IN phase separation occurs. This is in contrast to previous studies of mixtures of 
infinitely thin platelets and hard spheres jiol . \v\ . which reported isotropic-isotropic (JJ) 
demixing if the size ratio of the mixture is small. In our DFT treatment, we do not find // 
coexistence in the range of size ratios 0.2 < Rd/Rs < 10. 



B. Mixtures of colloidal platelets and ideal polymers 

We turn to the analysis of mixtures of colloidal platelets and ideal model polymers. Bates 
and Frenkel studied such mixtures using Monte Carlo simulations and perturbation theory 



26|. Following the ideas of Asakura, Oosawa and Vrij 27|, |28|, the polymers are modelled as 
freely interpenetrable spheres that cannot overlap with the platelets due to hard repulsion 
between the unlike species. A DFT for such mixture is obtained by linearizing the excess 
free energy, Eq. ([2D, around = (see Sec. Ill Bl for details). 

The bulk fluid phase behaviour of mixtures of colloidal platelets and non-adsorbing poly- 
mers is depicted in Fig. [3] for three different values of the size ratio Ro/Rp- The upper 
row shows the phase diagram in the pressure-composition plane. The polymer fugacity, 
Zp = exp{(]fip)/Vp, with fip the chemical potential of the polymers, as a function of the 
density of platelets is represented in the lower row. Here the chemical potential of species i is 
obtained from dF/dNi = fii. Note that the fugacity is independent of the (hence irrelevant) 
thermal volume Vp because the chemical potential, fip is shifted according to the particular 
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FIG. 3: Phase diagrams of mixtures of colloidal platelets and ideal model polymers for size ratios 
Rd/Rp = 2 (left), Rd/Rp = 5 (middle), and R^/Rp = 10 (right). Upper row: pressure, /SpR^, 
as a function of the composition of polymers, x. Bottom row: fugacity of the polymers, zpR^, 
vs density of platelets, ppiR^. Empty circles present critical points. Dotted lines represent triple 



points. Full circles in panel (bl) indicate Monte Carlo simulation data 



26|. 



choice of Vp. 

The main difference to platelet-colloid mixtures is that mixtures of platelets and ideal 
polymers display prominent demixing between two isotropic states when the size ratio is low 
enough. In panel (al) of Fig. |3]we show the results for R^/Rp = 2. There is strong demixing 
between two isotropic phases bounded by a lower critical point. The II phase separation 
occurs almost entirely at pressures below the IN pressure of the pure fluid of platelets. At 
higher pressures the II binodal merges with the IN binodal via a nematic-isotropic-isotropic 
triple point. This phase diagram differs significantly from that of the hard core mixture with 
the same size ratio shown in panel (al) of Fig. |2l The difference between both systems lies in 
the sphere-sphere excluded volume interactions. In a mixture of platelets and ideal polymers 
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the gain in accessible volume in the demixed // state over-compensates the loss of mixing 
entropy. This is not the case in mixtures of colloidal platelets and spheres for the size ratios 
considered. It is interesting to note that II demixing has been found by computer simulation 



4l|. 



in mixtures of platelets with finite thickness and non-adsorbing polymers 

In panel (bl) of Fig. [3] we show a comparison with Monte Carlo simulation data (full 
circles) 26| for the same mixture {Rd/Rp = 2). The agreement is reasonably good, although 
DFT underestimates the IN coexisting densities at low polymer concentration and also the 
fugacity of the polymers along both branches of the II binodal and hence the location of 
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the critical point. This trend is similar to what was found for the AOV model 

The middle (right) column of Fig. [3] shows the phase diagrams of a mixture with size 
ratio Rpi/Rp = 5 {Ro/Rp = 10). The behaviour is very similar to the corresponding 
mixtures of colloidal spheres and platelets, depicted in panels (a2) and (a3) of Fig. [2l Such 
similarity might be expected because by increasing the relative size of the platelets (in hard- 
core platelet-sphere mixtures) the platelet-platelet and platelet-sphere interactions dominate 
over the sphere-sphere interactions. Recall that we do not consider the freezing transition of 
hard spheres, which would dramatically change the phase behaviour between ideal polymers 
and colloidal spheres at densities above the freezing transition. 

Our results are in qualitative agreement with the experimental observations. For exam- 



ple, // coexistence has been found in suspensions of sterically stabilized colloidal 
platelets and non-adsorbing polymers (with R^/Rp ~ 2.9) by Lekkerkerker et al. 



_gibbsite 
The 



same study also showed a strong broadening of the IN two-phase region by increasing the 



polymer concentration. Nematic-nematic demixing has been recently found 



44| in a related 



system of mixtures of non-adsorbing polymers and positively charged platelets. 



C. FMT vs Parsons-Lee theory 

The geometry-based DFT of Sec. Ill Al contains contributions to the excess free energy 
that are of third order in density. In this section we compare the results from this theory to 
those from an extended Onsager theory with Parsons-Lee rescaling that explicitly includes 
only two-body interactions. For details about Parsons-Lee theory, see Sec. Ill Dl The phase 
diagram of a hard core mixture with size ratio Rpi/Rs = 2 according to Parsons-Lee theory 
is depicted in Fig. HI It is to be compared with panels (al) and (bl) of Fig. O where 
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FIG. 4: Bulk phase diagrams of a colloidal mixture of platelets and spheres according to Parsons- 
Lee theory, (a) Scaled pressure PpRjj vs composition of spheres x. (b) Density of platelets PdR% 
as a function of the packing fraction of spheres psvs- The size ratio is Rjj/Rs = 2. 

corresponding results obtain from FMT are shown. An immediately noticeable difference is 
in tlie location of the IN phase transition in the pure fluid of platelets. In this limit (x = 0) 
Parsons-Lee theory gives an unchanged result as compared to Onsager theory due to the 
vanishing packing fraction in the pure fluid of platelets. Hence Parsons-Lee theory clearly 
overestimates the coexisting densities (p^i?|, = 0.67 and p^R% = 0.86). It also predicts 
a much too high nematic order parameter at the IN coexistence, = 0.80 (not shown), 
compared to FMT and MC simulation. 

For Ro/Rs = 2 the topology of the phase diagram predicted by both theories is the same. 
According to Parsons-Lee theory the partitioning between isotropic and nematic phases is 
higher than in FMT. A similar trend was observed comparing the current FMT and Onsager 
theory in binary mixtures of hard platelets with different radius 12j. Parsons-Lee theory 
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predicts demixing between two nematic phases in mixtures with high size ratio, as does FMT. 
However, the minimal size ratio for phase segregation to occur between two nematics (using 
Parsons-Lee theory) is R^/Rs ^ 15, considerably higher than the threshold for FMT, which 
is R£)/Rs ~ 10. One important difference is that Parsons-Lee theory predicts // demixing 
for low size ratios Rd/ Rs ^ 1. Nevertheless, the II phase boundary is located at very high 
pressures and densities and is most likely metastable with respect to segregation between 
an isotropic phase rich in platelets and a solid phase rich in spheres. 



IV. CONCLUSIONS 



' mixtures of colloidal platelets and 
^ . The size ratio Ro/ Rs is the key 



We have investigated the fluid bulk phase behaviour o: 
spheres using a geometry-based density functional theory 
parameter that controls the behaviour of the system. In mixtures with 0.2 < R^i/Rs ^ 10 
we find only IN phase separation. We have identified two different mechanisms behind the 
IN phase transition. If the composition of spheres is low, the IN transition is driven by 
the excluded volume interaction between the platelets, similar to the mechanism in a pure 
fluid of platelets. In this regime, the transition takes place between two phases with low 
partitioning. As the composition of spheres is increased, the interaction between dissimilar 
species becomes dominant, and the (unfavourable) excluded volume interactions between 
spheres and platelets drives strong demixing between an isotropic sphere-rich phase and a 
nematic platelet-rich phase. 

Mixtures with high size ratio, R^/Rs ~ 10, also display demixing between two nematics 
at different compositions. However reliably describing such highly asymmetric mixtures is 
a very difficult for any theoretical treatment, see the comparison of FMT and simulation 
results for asymmetric binary hard sphere mixtures 



45 



46|. 



The possible existence of demixing between two isotropic fluids in mixtures of colloidal 
platelets and spheres remains an open question. The FMT approach does not predict II 
demixing. Previous studies 16|, ll7| reported II phase separation in mixtures were the size of 



both species is similar. Nevertheless, in both studies the authors neglected excluded volume 
interactions between platelets. This approximation is valid only if the platelets density is 
very low. In contrast our theory includes platelet-platelet interactions, modelled via up to 
cubic contributions in platelet density to the free energy. The bulk free energy for isotropic 
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47| . which 



states from the present DFT is the same as that from scaled-particle theory 
gives further confidence into our results. Carrying out a careful simulation study in order to 
shed further light on the existence of II demixing in the hard core mixture would be very 
valuable. 

We have not considered spatially inhomogeneous phases in this work. The pure fluid of 
hard spheres undergoes a fluid-solid phase transitions. The freezing transition, according 



to Monte Carlo simulations 
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49|, takes place at packing fractions psVs = 0.49 for the 



liquid and 0.54 for the solid (values which are higher than the packing fractions considered 
in the present work). Hence, we expect our fluid phase diagrams to be stable with respect 
to freezing of the hard spheres. For low size ratios (smaller than those analysed here) we 
expect an isotropic-solid demixing region in the phase diagram. Columnar as well as solid 
phases become relevant in systems of platelets with non-vanishing thickness at sufficiently 
high densities. Hence, it is possible that nematic-nematic demixing in mixtures with high 
size ratio will be metastable with respect to columnar ordering. Recent experiments by 



Lekkerkerker et al. on mixtures of charged gibbsite platelets and silica spheres 
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have shown a large isotropic-columnar coexisting region. Columnar ordering has also been 



studied theoretically in mixtures of board-like platelets and spheres in a cubic lattice 



50|. 



Smectic phases are not very common in discotic liquid crystals, but they could appear in real 



systems with finite thickness 



5l| . Recent experiments with mixtures of colloidal gibbsite 



platelets and spheres 52| show that also glass formation can occur at sufficiently high sphere 
concentrations. 

As our theory is mean-field in character, we expect the location of nematic-nematic 
critical point to change upon including neglected fiuctuations, say in computer simulations. 
Such (future) work could ascertain the stability of nematic-nematic demixing. Furthermore 
it would be interesting to consider the stability of the nematic-nematic transition upon 
altering the platelet-polymer interactions, in order to take into account penetrability of 
platelets and polymers. 

The results presented in this work can form the basis for studies of inhomogeneous 
platelet-sphere mixtures using DFT, such as the analysis of the effects of gravity |53|, as 
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well as wetting and confinement effects. 
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